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We study the mechanism how critical end points of first-order valence transitions are controlled by 
a magnetic field. We show that the critical temperature is suppressed to be a quantum critical point 
(QCP) by a magnetic field and unexpectedly the QCP exhibits nonmonotonic field dependence in 
the ground-state phase diagram, giving rise to emergence of metamagnetism even in the intermediate 
valence-crossover regime. The driving force of the field-induced QCP is clarified to be cooperative 
phenomena of Zeeman effect and Kondo effect, which create a distinct energy scale from the Kondo 
temperature. This mechanism explains peculiar magnetic response in Celrlns and metamagnetic 
transition in YbXCu4 for X=In as well as sharp contrast between X=Ag and Cd. 
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Quantum critical phenomena is one of the most inter- 
esting issues in condensed matter physics. Quantum crit- 
ical points (QCP) emerging when continuous-transition 
temperature of symmetry breakings such as magnetic or- 
ders is suppressed to absolute zero have been studied 
intensively [1, 2, 3], since the critical fluctuations induce 
unusual normal-state behaviors and even trigger other in- 
stabilities such as unconventional superconductivity [J- 

Valence transition is an isostructural phase transition, 
which is known to occur as 7-a transition in Ce metal @ 
and also in YbInCu4 [bI] characterized by a jump of the 
valence of the Ce and Yb ion. At the critical end point 
of the first-order valence transition (FOVT), the valence 
fluctuation diverges as diverging density fluctuations in 
the liquid-gas transition. When the critical temperature 
is suppressed by tuning material parameters and enters 
the Fermi-degeneracy regime, diverging valence fluctua- 
tions are considered to be coupled with the Fermi-surface 
instability. This multiple instability seems to be a key 
mechanism which dominates the low-temperature prop- 
erties of the materials includin g th e valence-fluctuating 
ions such as Ce and Yb 0, i, UU^- 

Actually, a remarkable increase of the superconducting 
transition temperature far from antiferromagnetic QCP 
in CeCuzXz (X=Ge, Si) [H [H d and Celrlng [H 
as well as linear-temperature dependence of resistivity 
observed in variety of Ce and Yb compounds [H, [l3|, 
[Hi [Hj seems to a sign of underlying influence 
of the proximity of the QCP of the valence transition 
(VQCP) [13. 

So far, to reduce the critical temperature of the va- 
lence transition toward absolute zero, intensive efforts 
have been made by chemical substitutions and applying 
pressure [l^. Usually, magnetic field also offers one of 
the efficient control parameters. However, it is highly 
nontrivial how the critical end point as well as the QCP 
is controlled by applying magnetic field, since valence in- 



stabilities are essentially ascribed to charge degrees of 
freedom, i.e., relative change of the 4f- and conduction- 
electron charges. 

In this Letter, we show that magnetic field is an ef- 
ficient parameter to change the critical end point of 
the FOVT to the QCP. Unexpectedly, we discover non- 
monotonic field dependence of the VQCP in the ground- 
state phase diagram, whose mechanism is clarified to be 
cooperative phenomena by Zeeman effect and Kondo ef- 
fect. We show that a metamagnetic jump in the magne- 
tization is caused by the proximity effect of the VQCP by 
magnetic field even in the intermediate valence-crossover 
regime. This mechanism explains a peculiar magnetic 
response observed in Celrlns [l^ as well as sharp 
contrast between YbAgCu4 and YbCdCu4 (l5[ • Our re- 
sults indicate significance of a distinct energy scale, which 
is characterized by the closeness to the VQCP as a key 
parameter for the valence-fluctuating materials. 

First we demonstrate how the critical end point 
is suppressed to absolute zero. Let us consider the 
Claudius-Clapeyron relation for the FOVT temperature 
Tv! ST^/Sh = — (mx — tomv)/(5'k — Smy), where m and 
S denote the magnetization and the entropy, respectively. 
Here, K represents the Kondo regime where the f-electron 
(hole) density per site (nf ) is close to 1 in the Ce (Yb) sys- 
tem, i.e., Ce^+(4f^) and Yb^+(4f^^), and MV represents 
the mixed- valence regime with (nf) < 1 [21j. Since the 
magnetization as well as the entropy in the Kondo regime 
are larger than those in the MV regime, as observed in 
the specific heat and the uniform susceptibility, it turns 
out that Tv is suppressed by applying h (see Fig. [TJa)). 
Then, the critical end point is suppressed to T = by ft.. 

Furthermore, the fleld dependence of in the zcro- 
tcmperature limit is also derived by using the above re- 
lation: For r ^ 0, the entropy shows the T-linear behav- 
ior in both the Kondo and MV regimes and Sk ~ Smv 
is approximated to be proportional to in case that 
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Tv is smaller than the characteristic energy scales in the 
Kondo and MV regimes. Noting the temperature inde- 
pendent mx — tomv, we have ST^/Sh — — Ci/Tv, leading 
to Tv = C2^/hv — h with constants Ci and C2, which 
explains well the observed behavior in the Ce metal |22] 
and YbInCu4 |23j (see Fig. [Ub)). We stress here that 
our analysis not only provides firm ground for small- 
Tv behavior by considering the coherence of electrons 
which is essential for low temperature, but also interpo- 
lates high Tv satisfying the relation [h/h^f + {T/T^f = 
1 [23, [13, [131 to zero temperature, since this relation 
was derived by assuming the isolated atomic entropy [2^ , 
which is justified only at the high temperature regime. 
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FIG. 1: (color online) (a) Schematic phase diagram, showing 
the FOVT surface in the T-P-h space, where P represents 
a control parameter (e.g. pressure, chemical concentration, 
etc.). The critical end point(CEP)s form a continuous tran- 
sition line touched at T = as the QCP. (b) FOVT line 
{h/Kf + {T/T^f = 1 [ll,!!!,!!! in the T-h plane for a fixed 
P, corresponding to the dashed line in (a) . 



Although we have shown that the h dependence of the 
critical end point can be understood from the viewpoint 
of the free-energy gain by the larger entropy in the Kondo 
regime, it is highly nontrivial how the VQCP is controlled 
by /i at T = 0. To proceed our analysis, we introduce a 
minimal model which describes the essential part of Ce 
and Yb systems in the standard notation: 



H — He + Hf + ffhyb + -fft/fc 



where Hr 



(Si 



Sf 



(1) 



^c/fc — ^fc 12i^=i^i^i- The Ufc term is the Coulomb 
repulsion between the f and conduction electrons, which 
is considered to play an important role in the valence 
transition [3, [13, [IE [13, [13i [13 ■ For example, in the case 
of Ce metal which exhibits the 7-a transition, the 4f- 
and 5d-electron bands are located at the Fermi level [1^ . 
Since both the orbitals are located on the same Ce site, 
this term cannot be neglected. In the case of YbInCu4, 
Hu^^ also plays a crucial role for the FOVT in the hole 
picture of eq. ^ [28] . 

We apply the slave-boson-mean-field theory [3] to 
eq ([T]). To describe the state for U — 00, we 
consider V fj^hiCicr instead of V fj^Cia in eq. ([1]) by 
introducing the slave-boson operator bi at the i-th 
site to describe the f*^ state and require the con- 
straint (So-^L + - 1) with Xi being the 
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Lagrange multiplier. For Hu^^ in 
ploy the mean-field decoupling as 
ncn\ — i^'^^iUc- By approximating mean fields as uni- 

set of 
and 
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form ones, i.e., h = {bi) and A = A^, the 
mean-field equations is obtained by d{H)/d\ = 
d{H)/db = 0: A = fiE^J-fiEL) 



\b\' = jwJ: 



^kcr,± 



1± 



f{E, 



kcr) 



and the 



following equation holds for the total electron num- 
ber: Uf + = Ek. [fiK.) + fiKa)] Here, 



f{E) is the Fermi distribution 
the lower (— ) and upper (-I-) 
the quasi particle with spin cr. 



function and E^^ are 
hybridized bands for 



kcr 



— £e 
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respectively: 

Sfo- + £Wa ± V {£fa - Ska)"^ + , whcrC Eko 

and V are defined by Eko- = £k + UfcUf 
Ef -f A -I- C/fcfic - and V = V%. The dispersion of 
the conduction electrons is taken as Ek = k^/(2TO) — D 
with D being the bottom of the conduction band and the 
density of states iVo(e) is set to satisfy the normalization 
condition, J_jj d£NQ{e) = 1 per spin in three dimension. 
We take D as the energy unit and show the results for 
V = 0.5 and the total filling n = (nf + nc)/2 = 7/8. 

At h — and T = 0, a jump in fif appears as a func- 
tion of £f for large Ufc, which indicates the first-order 
quantum phase transition, since a jump in Hf results 
in the level crossing of the ground states by the rela- 
tion fif = d{H)/def. The FOVT between the Kondo 
state and the MV state is caused by Ufc, since a large 
Ufc forces the electrons to pour into either the f level or 
the conduction band. The QCP in the Ef-Ufc plane is 
identified to be {sf^^, U^'^^) = (0.356, 1.464), at which 
the jump in fif disappears and the valence susceptibility 
Xv = —d^{H)/def = —dfif/dsf diverges. The charac- 
teristic energy scale of the system, the so-called Kondo 
temperature, which is defined as Tk = Sfa- — M is esti- 
mated to be Tk = 0.074 at the QCP. 

Note that the surface of the FOVT exists in the pa- 
rameter space of £f, Ufc and V and a trajectory line is 
drawn in the space for corresponding experimental pa- 
rameter such as pressure (i.e., P in Fig.l). It is also noted 
that the valence instability is considered to be coupled 
to the phonon degrees of freedom. Hence, in our model 
([1]) the effect of the hybridization also plays an impor- 
tant role, which might share common aspects with the 
Kondo- volume-collapse scenario for Ce metal 

A remarkable result is found under the magnetic field 
in the valence-crossover regime for Ufc < C^f^^^- Fig- 
urelUa) shows the magnetization m = J2i{^i^ + Si^)/N 
vs. h for {ef,Ufc)= (-0.354,1.458) (thin line) and 
(—0.349, 1.442) (bold line), indicating that the metamag- 
netism defined by the diverging magnetic susceptibility 
X = dm/dh — oo emerges at /i = ft.m = 0.01 and 0.02, 
respectively. To clarify the origin, we have determined 
the FOVT line as well as the QCP under the magnetic 
field. The result is shown in Fig. [3] It is found that 
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FIG. 2: (color online) m-h curve (a) for (sf, Uic) ~ 
(-0.354, 1.458) (thin hne) and (-0.349, 1.442) (bold line), and 
(b) for Ef ranging from —0.36 to —0.32 at C/fc = 1.42. In both 
cases, _D = 1, V = 0.5 at n = 7/8. 



the FOVT line extends to the MV regime and the lo- 
cation of the QCP is shifted to the smaller Ufc and |ef| 
direction, when h is applied. This \ow-h behavior of the 
FOVT line agrees with the low-temperature limit of Tv 
discussed above. 

We show the m-h curve at U[c = 1.42 for Sf ranging 
from -0.32 to -0.36 in Fig. ^h). Tk at ft. = is esti- 
mated as 0.0353, 0.0873, 0.1346, 0.1611, and 0.1823 for 
Ef = -0.36, -0.35, -0.34, -0.33 and -0.32, respectively 
(see crosses in Fig.|31). From these results, the mechanism 
is understood as follows: At ft, = 0, Tk is originally large 
for Ef = —0.32 and —0.33, since the system is in the MV 
regime. However, by applying ft, the QCP approaches, 
which makes Tk reduced, since the system is forced to be 
closer to the Kondo regime by ft. At the magnetic field 
ft — ftm where the QCP is reached, the metamagnetism 
occurs with a singularity 6m ^ 5h^/^ [soj as shown in 
Fig. mja). On the other hand, for Sf = —0.35 and —0.36, 
the QCP is not approached, so that no metamagnetism 
appears. 

An intriguing result is found in Fig.[3l which exhibits a 
non-monotonic ft dependence of the QCP: As ft increases, 
the QCP shows an upturn around ft = 0.04, which is 
comparable to Tk at the QCP for ft = 0. The upturn of 
the QCP has been also confirmed for a constant density 
of states No{e) = 1/(2D), suggesting that this behavior 
appears irrespective of details of the band structure. 

The non-monotonic behavior can be understood from 
the structure of the valence susceptibility Xv, which 
is given essentially by the RPA as discussed in ref.l7. 
Namely, it is given as Xv(g) ~ Xfc\<l)/[^ - C^fcXfcH?)], 
where x^^' is the bubble diagram composed of f and 
conduction electrons. In the Kondo regime (ft < Tk), 
where f electrons have a dominant spectral weight at 
around e ^ — et with width A ~ ttV'^N{sf), Xfc'' is 
estimated as Xfc'' ~ Vkd and is shown to be an in- 
creasing function of ft. Therefore, U^'^^ is decreas- 
ing as ft is applied until it reaches around ft ~ Tk, 
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is also decreasing. For ft > Tk, 
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FIG. 3: (color online) Ground-state phase diagram in the 
plane of [/& and et for D = 1, 1/ = 0.5 at n = 7/8. The FOVT 
line with a QCP for h = 0.00 (open triangle), 0.01 (filled 
triangle) 0.02 (inverse triangle), 0.03 (filled star), 0.04 (filled 
diamond), 0.05 (filled square) and 0.06 (open square). The 
shaded line connects the QCP under h, which is a guide for 
the eyes. The dashed line with crosses represents Ui^ = 1.42 
(see text). 



xlc^ ~ 1/A < l/|£f| with a help of shift of the f level 
towards the Fermi level, i.e., Sf ^ e{ + UfcS{nc) {S{nc) 
being the change of number of conduction electrons per 
site due to entering the MV regime), so that U^'~^^ is 
larger than that of C/f^^^(ft Tk). 



mass enhancement (~ 1/^) is quickly suppressed and 
the MV regime is approached. Then, x^^-* is given as 



We find that ftm corresponds to the difference of Tk at 
the QCP between ft = and ft 7^ 0, i.e., ft^ ~ AT^°^ = 
T]^^^(ft 7^ 0) - T^^^(ft = 0). This indicates emergence 
of a new energy scale distinct from Tk, which is char- 
acterized by the closeness to the VQCP. Furthermore, 
the proximity of the VQCP to the intermediate valence- 
crossover regime under ft yields emergence of the meta- 
magnetism as well as a jump in m even without showing 
the temperature-driven FOVT at ft = 0. 

This mechanism explains a peculiar behavior observed 
in Celrlns, which shows a jump in the m-h curve at 
42 T, but does not show the first-order transition in any 
physical quantities when temperature is changed [Tol . |20| | . 
Namely, this is naturally understood if Celrlns is located 
inside the enclosed area of the QCP line for ft 7^ in 
Fig. [3l This provides a key to resolve the outstanding 
puzzle about the origin of the superconductivity of this 
material, whose transition temperature increases even 
though antiferromagnetic spin fluctuation is suppressed 
under pressure [14]. Since superconductivity is shown to 
appear near the VQCP 0, [lo|, this view verifies a new 
scenario that the proximity of the VQCP is the origin of 
the superconductivity. 

Our result also explains why the metamagneic in- 
crease appears in the m-h curve in YbAgCu4, but not in 
YbCdCu4, in spite that both have the comparable TV's 
estimated from the T ^ magnetic susceptibility [15i |. 
The relation ftm ~ AT^*"^^ clearly indicates that different 
ft,n can appear for the same Tk according to the location 
in Fig. [3l Namely, the sharp contrast can be understood 
if YbAgCu4 and YbCdCu4 are located in the valence- 
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FIG. 4: (color online) Magnetization process for Sk ~ 
-2cos(fc), V = 1, U = 10*, Si = -1 and [7fc = 1 at n = 7/8: 
(a) m-h curve (filled circle). An arrow indicates the metam- 
agnetic transition. Inset: et dependence of nf extrapolated to 
the bulk limit for Ufc ~ (open circle), f/fc = 1 (filled trian- 
gle) and U{c ~ 2 (open triangle). An arrow indicates Sf — —1. 
Lattice structure used in the calculation, (b) fif (filled trian- 
gle) and He (filled square), nif (open triangle) and ttIc (open 
square). In (a) and (b) dashed lines represent h = hm- 

crossover regime for Ufc < U^'^^ and AT^*~^^ for the 
former is smaller than that for the latter. Furthermore, 
the evident FOVT as a function of T as well as a 
jump in the m-h curve [3l| observed in YbInCu4 is also 
understood if this material is located in the MV regime 
for C/fc > C/f^^^ in Fig. [3] 28]. This is consistent with 
recent band-structure calculations [3^]. 

To examine the mechanism more precisely, we have ap- 
phed the density-matrix-renormalization-group method 
to eq. ((D) in one dimension (ID). Since valence fluctu- 
ations are basically ascribed to the atomic origin, the 
fundamental properties are expected to be relevant even 
in ID 10]. We show here the results for Sk = — 2cos(fc), 
V = 1, U = 10"* at n = 7/8 on the lattice with = 40 
sites for conduction electrons illustrated in the inset of 
Fig- which may be regarded as a ID mimic of 



Celrlns and YbXCu4. For h = 0, fif vs. £f for U[c = 0.0, 
1.0 and 2.0 is shown in the inset of Fig. HJa). As Ufc 
increases, the change in nf becomes sharp. We show 
here the magnetization m in the MV state of which nf at 
/i = is indicated by an arrow in the inset of Fig. HJa). 
A plateau appears at m = 1 — n = 1/8, which is ex- 
pected to disappear if we take more realistic choice of 
parameters, e.g. the momentum dependence of V and 
Ef. The main result is that metamagnetism emerges as 
indicated by an arrow. The increase of nf simultaneous 
with the decrease of he a,t h — shown in Fig. Sfb) 
reveals that this is caused by the field-induced extension 
of the QCP to the MV regime. Namely, these results in- 
dicate that the mean-field conclusion is not altered even 
after taking account of quantum fluctuations and elec- 
tron correlations. 

To explore further the nature of this metamagnetism, 
we have calculated nif — J2ii'^i^) ^^'^ = 
J2i{^i^) and we find that TOc decreases at h — h^, 
while mf increases as shown in Fig.dfb). Since the Kondo 
cloud is still formed even at /i = h^, i.e., (Sf • Sf) < 0, 
a remarkable decrease of (Sf^) is ascribed to the field- 
induced Kondo effect, which is a consequence of the en- 
ergy benefit by both the Kondo effect and the Zeeman 
effect. Although this mechanism itself has been known 
to occur in the Kondo regime (33| . the present finding 
is that such a mechanism works in the MV regime as a 
driving force of the field-induced VQCP. 

In summary, we have clarified the novel mechanism 
of the magnetic-field dependence of critical points of va- 
lence transitions and we believe that the proximity of the 
VQCP revealed in this study provides an essential ingre- 
dient to understand the T-P-h phase diagram of the Ce 
and Yb systems. 
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